Ninit = [1.0; 0.0; 0.0]
c = [2.5; 2.5; 2.5]
R = 8.13
V0 = 1.0
T0 = 1000.0
Vreact = 1.0
Treact = 1000.0
s0 = [248.462; 243.507; 228.087]
